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ABSTRACT 

We present high dispersion (R ~ 16,000) optical (3900 — 8700 A) spectra of 
390 stars obtained with the Palomar 60-inch Telescope. The majority of stars 
observed are part of the Spitzer Legacy Science Program The Formation & Evolu- 
tion of Planetary Systems. Through detailed analysis we determine stellar prop- 
erties for this sample, including radial and rotational velocities, Li I A 6708 and 
Ha equivalent widths, the chromospheric activity index R'^ki ^'^d temperature- 
and gravity- sensitive line ratios. Several spectroscopic binaries are also identified. 
From our tabulations, we illustrate basic age- and rotation-related correlations 
among measured indices. One novel result is that Call chromospheric emis- 
sion appears to saturate at vsmi values above ~ 30 km s~^, similar to the well 
estabhshed saturation of x-rays that originate in the spatially separate coronal 
region. 

Subject headings: stars: fundamental parameters — stars: activity — stars: pre- 
main sequence 



1. Introduction 

We have obtained high dispersion spectra of several hundred solar type stars within 
20-160 pc of the Sun in order to determine fundamental stellar properties. Our work is in 
support of the Spitzer/ Legacy Program "Formation and Evolution of Planetary Systems" 
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(FEPS; iMeyer et al. II2006I ). which aims to understand the evolution of circumstellar dust 
from the primordial planet-building phase at an age of ~ 3 Myr to mature debris disk sys- 
tems up to an age of ~ 3 Gyr. While data from the Spitzer Space Telescope probe the 
dust characteristics, a wide range of ancillary ground-based observations are required to 
estimate stellar properties. High dispersion spectroscopy, in particular, is a valuable tool 
for determining radial and rotational velocities, effective temperatures and surface gravities, 
abundances, and chromospheric activity diagnostics. These in turn provide useful constraints 
on stellar masses and ages. The stellar parameters derived from these spectroscopic obser- 
vations thereby permit multiple axes of investigation for the infrared data from Spitzer, such 
as whether dust disk characteristics can be (anti-) correlated with stellar mass, signatures of 
youth, and/or stellar multiplicity. 

In Section 2 we describe the observed sample and the spectroscopic observations. In 
Sections 3 and 4 we summarize the spectroscopic reductions and the methods used to extract 
spectroscopic properties. In Section 5 we discuss the spectroscopic binaries identified in the 
sample and illustrate basic age related correlations among the measured properties; the 
relation and apparent saturation of chromospheric Call emission versus projected rotational 
velocity is highlighted. The analysis is meant to be representative rather than complete. The 
spectroscopic properties are presented primarily to assist more detailed studies of individual 
FEPS stars. 



2. Sample and Observations 

Our primary sample for high dispersion spectroscopic observations consists of actual 
FEPS targets as well as many candidate targets eventually dropped from the FEPS pro- 
gram. The source list for FEPS is comprised of young near-solar analogs that range in 
mass from 0.7-1.5 M0 and span ages between 3 Myr and 3 Gyr. This source list was drawn 
from three samples of young solar-analogs. The first sample was assembled from the D. R. 
Soderblom (2000, private communication) volume-limited (< 50 pc) spectroscopic study of 
stars with B-V colors between 0.52 and 0.81 mag (spectral types F8-K0) and Mv mag- 
nitudes within 1.0 mag of the solar- metallicity zero- age main sequence; approximate ages 
are provided from the chromospheric activity index R'hk- Since this nearby field sample 
primarily consists of stars older than a few x 10^ yr, a second sample was assembled with an 
targeted age range of 3-300 Myr. Stars in this sample were identified based on observational 
investigations to identify young main-sequence and pre-main sequence solar-type stars with 
strong x-ray emis sion, high lithiuna abu ndances, and kinematics appropriate for the young 



galactic disk (e.g. iMamajek et al. II2002I ). The third sample consists of stars in nearby well- 
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studied clusters, including 4 open clusters [IC 2602 (55 Myr), Alpha Per (90 Myr), Pleiades 
(125 Myr), Hyades (650 Myr)] and 3 younger T Tauri age clusters [R Corona Australis (3 
Myr), Upper Scorpius (5 Myr), Upper Centaurus Lupus (10 Myr)]. Stars in these clusters 
serve to benchmark our field star selection criteria by providing sub-samples nearly identical 
in age, composition, and birth environment. These three parent samples and the procedure 
for downselec t ing to the final FEPS sample observed with Spitzer are described more fully in 
Meyer et al. I (120061 ). In addition to stars in this primary sample, a handful of other young 
and/or chromospherically active stars (e.g. RS C Vn) were observed fo r comparison, and 
~ 20 stars with precisely known radial velocities (INidever et al. 1 12002| ) were observed for 
calibration and analysis purposes. 

High dispersion spectroscopic observations were obtained for 363 stars among the sample 
described above; 284 of these stars are among the 326 stars that comprise the final FEPS 
program as executed with Spitzer (the vast majority of FEPS stars not observed are in the 
southern hemisphere and thus inaccessible from Palomar Observatory). Fifty-two of the 
363 stars observed are members of the 3 open clusters and 31 stars are members of the 
3 young clusters. For calibration purposes, we also observed 13 giant stars with spectral 
types ranging from G5 to M2, 5 early- to mid-M dwarf stars, and 9 early A type stars with 
nearly featureless continua. Tables 10 and 2 list the observed sample (excluding non-solar 
type calibration stars), ordered by right ascension; Table 1 contains spectroscopic single 
stars and candidate single-lined spectroscopic binaries (for which stellar properties could be 
determined) while Table 2 contains known or newly discovered double-lined spectroscopic 
binaries. Coordina tes for all stars are from t he 2MASS All-Sky Cata log of Point Sources 
(ICutri et al. II2003I ). as interpreted by VizieR (lOchsenbein et al. Il2000l ). 



Spectroscopic observations were obtained on 9 observing runs betwe en July 2001 and 
June 2003 with the Palomar 60- inch telescope and facility spectrograph (IMcCarthy Ill988l ) 
in its echelle mode. This instrument is no longer available. A 1'.'43 x 7'.'36 slit was used, 
yielding spectra spanning from approximately 3900 A to 8700 A, with some small gaps 
between the redder orders. While the 2-pixel resolving power of this instrument is ~ 19, 000, 
the achieved resolving power was typically ~ 16, 000 because of image quality issues. The 
wavelength coverage was chosen so as to include both the Ca II H&K features at the blue end 
and the Call infrared-triplet at the red end. The spectrograph images the echelle spectra 
onto an 800x800 CCD; the consequence of projecting the broad wavelength coverage onto 
this modest sized CCD is that the orders are closely spaced, especially at the red end. The 
detector has a gain of 1.5 and read noise ~ 12 electrons. During each night of observations. 



^The complete version of Table 1 is available online. 
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series of bias exposures and flat field exposures were obtained to effect standard image 
processing. Because of the large spectral range of the instrument, the flat field images 
were taken with two separate incandescent lamps, one appropriate for illumination of the 
blueward range of the detector (at 60 seconds exposure time) and the other appropriate 
for the redward range (5 seconds exposure time). For wavelength calibration, spectra of 
Thorium-Argon lamps were obtained (45 seconds exposure time) generally before or after 
each on-sky group of target exposures. 

Each target observation consisted of between 1 and 3 consecutive exposures. The dates 
for these observations are listed in Tables 1 and 2; multiple epoch observations are listed 



separately. At least 3 stars with precisely known radial velocities from lNidever et al. I ((2002) 
were observed per observing run. 

One feature of this instrument is that the optimal focus for the bluest orders is not 
optimal for the reddest orders, so subjective compromises had to be made in focusing the 
spectrograph. Additionally, the spectroscopic focus was temperature sensitive, and therefore 
had to be tuned periodically throughout the night to achieve the best spectral resolution. On 
the night of 2002 Feb 2, many spectra were obtained at an exceptionally poor focus setting. 
When poorly focused, the spectrograph produces spectra that are stigmatic - spectral lines 
skewed relative to the dispersion direction which degrades resolution (by an amount that is 
wavelength dependent). To help assess how poor focus may compromise the analysis (e.g. 
vsini measurements), on the night of 2002 Sep 18 we observed the radial velocity standard 
HD 164992 over a broad range of focus values. 



3. Spectroscopic Reductions 

3.1. Image Processing 

Raw data frames were processed using tasks in IRAF and custom IDL scripts. First, 
the detector bias was removed. Although the bias pattern appeared to be stable during 
the course of an observing run, the absolute bias level drifted substantially on time scales 
of hours, changing by as much as 30% over the course of a night. Bias subtraction was 
performed, therefore, in a two-step process. First, the bias sequence obtained on each night 
was median combined using the IRAF "zerocombine" task, rejecting the highest and lowest 
valued bias images. Second, since the detector bias level for a given image may have drifted 
from the median value determined at the beginning or end of the night, a bias scaling factor 
was also determined. This was done by measuring the median pixel value within a specified 
area in the most poorly illuminated, inefficient portion of the array, located well below 3600 
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A where useful spectra could not be extracted; even the brightest, bluest stars observed on 
our program showed essentially zero counts above bias in this portion of the array. For each 
target image, the median bias image from the beginning or end of the night was scaled by 
this factor and then subtracted. 

Next, images were trimmed to exclude two unilluminated columns at the rightmost edge 
of the array. Cosmic ray events were removed using IRAF's "cosmicrays" task and fairly 
liberal parameter settings that restricted the corrections to only bona fide cosmic rays. 

The blue and red flat field sequences from each night were median combined using 
IRAF's "flatcombine" procedure with sigma clipping. The blue and red medians were then 
stitched together approximately 1/3 of the way up the detector (at line 250) which has 
the effect of distorting the one (curved) spectral order which intersects this line. IRAF's 
"apfiatten" task was used to normalize the combined flat for each night by fltting its intensity 
along the dispersion using a low order (typically 3) polynomial while setting all pixels outside 
an order to unity. Flat fielding was performed by dividing each image by the normalized fiat 
field image using IRAF's "ccdproc" task. 

Cosmetically, the images contain "tadpole streaks" of high count levels, comprised of 
short "heads" along the dispersion direction and then decaying count levels along long, mostly 
single pixel, "tails" perpendicular to the dispersion direction which thus affect multiple 
orders. The number of such defects per image was related to exposure time. A custom 
procedure was developed to crawl down each column comparing values on either side of the 
column and thus identifying these features. The pixels in the tail of the streak and a small 
square around the tadpole head were interpolated over. 

The close spacing of the spectral orders did not permit measurement or removal of 
scattered light, an effect common in echelle spectrographs. However, comparisons of our 
extracted spectra with spectra previously obtained with other instruments (Section 3.2) and 
of measured equivalent widths with those from the literature (Section 4.1) suggest that the 
majority of excess fiux caused by scattered light is removed during background subtraction. 

3.2. Spectral Extraction and Wavelength Calibration 

One dimensional spectra were extracted from the processed two dimensional images 
using IRAF's "doecslit" task. Sixty orders were identified over the 800 x 800 pixeP array. 
A moderate (9th) order trace was employed and the spectra in each order were summed 
over 3-5 pixels perpendicular to the dispersion. Variance weighting was used for "optimal" 
extraction. The small separation between orders made removal of background and scattered 
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light especially challenging. After exploring a variety of methods for accomplishing this, the 
most robust method identified was simply to subtract the average minimum value within 
2 background regions, one on each side of the spectrum. Aperture definition, background 
region definition, and order tracing were attended as opposed to automated processes. 

To test the accuracy of the adopted background/scattered hght subtraction, Jeff Valenti 
kindly provide d a careful compa rison of our spectra with spectra obtained with the HIRES 



spectrograph (jVogt et al. Ill994l ) on the W. M. Keck I Telescope. The widely spaced orders 
of HIRES permit accurate background and scattered light corrections. These comparisons 
demonstrated that the average scattered and/or background light remaining in our spectra is 
approximately 4% at 3900 A, and drops roughly linearly to 0% at 6200 A. The implication is 
that equivalent widths of features at blue wavelengths may be diminished by a small amount. 
Rather than attempt to correct for this in the reduction process, we account for this in the 
relative calibration of properties extracted from blue wavelength features (Sections 4.3 and 
4.4). 

Although multiple Thorium-Argon spectra per night were obtained (often one for each 
star observed), in practice we established wavelength calibration as follows. The extraction 
trace defined by a single bright star on each night was applied to its corresponding Thorium- 
Argon exposure to extract a one-dimensional wavelength reference spectrum. A wavelength 
solution was determined using the "ecidentify" task in IRAF, and this solution was applied 
to all spectra for a given night. This established wavelength solution was then interpolated 
onto a log-linear scale. 

From analysis of spectra taken on the same night, small spectral shifts of ~ 0.1 A were 
common, as inferred from comparisons of the location of telluric absorption lines. These shifts 
were hkely due to flexure of the instrumental optics with telescope pointing, and consequently 
resulted in wavelength solution offsets of this order for many spectra. As discussed below, 
this wavelength offset is accounted for in determining radial velocities by cross-correlating 
telluric features. 



4. Spectroscopic Properties 

Figure [T] shows portions of the echelle spectra for 7 stars near the Call H&K, Li I A6708 
A, and Ha features. To help assess the quality of the spectra, two signal-to-noise ratios 
(SNRs) are estimated, one at ~ 6700 A and one at ~ 3950 A (labeled as Li and HK, 
respectively). The SNR at 6700 A is determined from the gain-corrected median signal 
within the central portion of this order. The SNR at 3950 A is determined similarly, but is 
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based the average of the values determined for the two adjacent orders containing the Ca II 
H (3968 A) and K (3934 A) features, and thus represent an approximate average over the 
continuum and broad absorption. These estimates are vahd above values of > 20, in which 
the noise is dominated by photon statistics. Lower values may actually over-estimate the 
SNR, because of the increased relative error in the background, scattered light and bias level 
subtractions (Section 3.1). The SNR distribution at 6700 A is shown in Figure [2l 

The following sections described the procedures used to extract stellar properties from 
these spectra. In comparisons with previous measurements, values determined in this study 
are referred to as P60 values. Specifically, equivalent widths (EWs) of Li I A6708 A and Ha 
are measured in Section 4.1, radial and rotational velocities are determined in Section 4.2, 
temperature and surface gravity indices are discussed in Section 4.3, and measured S values 
and calculated R'hx indices are presented in Section 4.4. 



4.1. Equivalent Widths of Li I 6708 A and Ha 

Equivalent widths of the Li I A6708 A doublet (A6707.76, A6707.91), which is spectro- 
scopically unresolved in our observations, and Ha (A6562.8 A) are measured for all spectra 
using the splot task within I RAF. For measurement of the Li I absorption feature which is 
located near the peak of the blaze function in these spectra, the local continuum is deter- 
mined by fitting a legendre polonomial to a ~ 30 A portion of the spectrum, excluding values 
less than 2a below the continuum (e.g. stellar absorption features) or more than 4a above 
the continuum (e.g. noise spikes) in this fit. The order of the polynomial is varied (typically 
from 4 to 7) to best match the local continuum. With the continuum defined in this way, the 
EW values are determined by fitting a Gaussian profile to the absorption feature. Multiple 
measurements on a single night are averaged, while separate epochs are listed independently. 
If no feature is discernable in the spectrum, upper limits are determined from the noise in 
the local continuum. Table 1 lists the measured values. 

The measured Li I EW values include a contribution from an Fel line at 6707.441 A. 



Soderblom et al. I (119931 ) find that the strength of this temperature sensitive Fel feature 
varies with B-V color as EW[Fe I A6707.441] = 20(5 -V)-3 mA. Since the S - ^ values 
of most stars are known or can be estimated (Section 4.4), the contribution of this Fel line 
to the Li I EW values can be removed. 

To assess possible systematic uncertainties caused by the subjective aspects involved in 
measuring EW values (e.g. continuum determination) or instrumental effects, we compare 
our measurements to previous values determined from high dispersion spectroscopy. Of the 
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Fig. 1. — Portions of the extracted spectra. The top 3 spectra are G5 stars with decreasing 
levels of Call H&K emission (logi?^^ = —4.13, —4.59, —5.01). The next 3 spectra are 
Pleiads (and thus coeval), with spectral types of F5 (HII 1200), G2 (HII 250), and KO (HII 
2106). With a SNR of only ~ 18 in the Call region, HII 2106 is typical of the poorest 
quality spectra from which log R'hk values can be extracted; strong core emission can be 
seen in this case. The bottom spectrum shows the spectroscopic binary HD 78899, whose 
features are clearly doubled. 
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Fig. 2. — Signal-to- Noise distribution of all spectra obtained (includes multiple measure- 
ments and stars observed for calibration) of the 390 unique targets. 
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spectr oscopically single stars observed here, 21 Ple iads were observe d by ISoderblom et al 



(119931, R = 50,000), 19 Hyads were obs erved bv iThorburn et al I fll993l . R = 26,000 



32,000), 26 field stars were observed bv IWichmann et al. I (120031 . R = 40,000), and 28 



field stars were observed by IStrassmeier et al. I (l2000l . R = 25,000). The uncertainties in 



the EW[Li]s from these previous studies are typically 5-10 mA. Of these 94 stars, 19 were 
observed twice by us, yielding a total of 113 measurements for comparison. Figure [3] shows 
the difference between the P60 and previous EW values; all values are as measured and 
thus uncorrected for the contamination of the Fe I line. Overall the agreement is good. The 
average and/or median difference for any of the 4 studies is ±15 mA, with standard deviation 
of differences that range from 16 - 21 mA. The average difference with the entire comparison 
sample is -8 mA with a standard deviation of 20 mA. We conclude that any systematic bias 
in our measurements is < 10 mA. We adopt a uniform uncertainty of 0.02 A. 

For Ha EWs, considerable care had to be taken both in normalizing to the local con- 
tinuum and in fitting a functional form to the line profiles because of the broad wings. The 
local continuum was determined by fitting a legendre polonomial generally of order 10 to 
the entire spectral order, excluding values less than 2a below or more than 4ct above the 
continuum. The procedure varied slightly as dictated by the noise level in the data, by the 
line breadth (because of, for example, rapid rotation), or if Ha was found to be in emission. 
Both the sigma-rejection limits and the fitting order were modified in these cases to opti- 
mize the continuum fit. The EW values were determined by fitting a Voigt profile to the 
absorption feature; in a few cases of dramatically filled in Ha, direct integration was used to 
determine the EW rather than function fitting. EWs are listed in Table 1; negative values 
indicate Ha emission. 

Four stars (HD 143006, RX J1842.9-3532, RX J1852.3-3700, LH98 196) show especially 
strong Ha emission (> 10 A); their Ha profiles are shown in Figure HI In these cases, the 
strength and breadth of these emission-line profiles are greater than that expected from 
chromospheric activity alon e, and more consistent with that expected fro m the accretion of 
circumstellar material (e.g. iMuzeroUe et al. Ill998l : IWhite fc Basri II2003I ). We suggest that 
these stars are accreting. 

For the 8 stars identified as double-lined spectroscopic binaries, EW[LiI]s and EW[Ha]s 
for the individual components were determined by fitting two Gaussian profiles to the doubled 
features, when possible. For these pairs, we also measured EWs of the temperature sensitive 
Cal A6717 A feature, which can be used to identify better the primary and the secondary. 
The EW values for the binaries are listed in Table 2. 
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Fig. 4. — Ha profiles for the 4 strong emission-line stars normalized by the local continuum; 
the Ha equivalent width for each star is given in parentheses. A velocity width of 300 km 
s~^ is indicated. 
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4.2. Radial and Rotational Velocities 

Radial velocities and projected rotational vel o cities are determined via a cross-correlation 



analysis similar to that of IWhite fc Hillenbrand I (120041 ) . The analysis is restricted to the 36 



orders spanning from ~ 4000-7500 A, excluding the 3 orders containing the telluric B-band, 
Ha, and Nail D features; orders at longer wavelengths have poorer wavelength solutions and 
those at shorter wavelengths typically have low SNRs. From this restricted range, orders 
having an average SNR greater than ~ 25 are then cross-correlated with the spectral orders 
of at least 3 sl owly rotating compar ison standards that have radial velocities accurate to 0.3 



0.4 km s ^ (INidever et al. 1120021 ). The radial velocity, relative to each standard, is then 



determined from the average velocity offset measured from all orders. 

To determine the heliocentric radial velocities of the stars, the relative radial velocity 
must be corrected for barycentric motions and possible errors in the wavelength solution of 
either the standard or target spectrum. The former is determined using the rvcorrect task in 
IRAF. The latter correction is determined by cross-correlating the telluric A- and B-bands 
of the standard and the target star; the telluric correlations would yield an offset of zero for 
perfect wavelength calibration, but are typically 1-3 km s~^ for the spectra analyzed here. 
The multiple heliocentric radial velocity estimates, determined from multiple radial velocity 
standards, are then averaged to find the final radial velocity. 

Uncertainties in the radial velocities are estimated from a combination of statistical 
and empirical error estimates. First, the uncertainty in the radial velocity relative to each 
stardard is assumed to be the uncertainty in the mean value of all orders used, combined 
with the uncertainty in the radial velocity of the standard (typically 0.3 km s~^). The 
radial velocity estimates from all standards are then combined using a weighted average; the 
resulting statistical uncertainties are typically small (0.1 - 0.3 km s"^). As an empirical check 
on these uncertainty estimates, the radial velocity of each standard is determined from all 
other radial velocity standards observed during that run and compared to its assumed value. 
The average agreement of the radial velocity standards is typically a couple times larger 
than the statistical uncertainties, except in low SNR cases where the statistical uncertainty 
dominates. This suggests there may be systematic effects unaccounted for in the statistical 
uncertainty estimates. Thus, we combine the statistical uncertainty with this empirical 
uncertainty to determine final (conservative) radial velocity uncertainty estimates given in 
Table 1. Figure E] shows the distributions of radial velocities for the field, open cluster, and 
young cluster samples. 

Rotational velocities (wsim) for the observed targets are determined from the width of 
the peak in their cross-correlation with a slowly rotating standard. Specifically, the peak 
is fit by a parabola with the fitted width at half-peak maximum. To convert this width 
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Fig. 5. — Distributions of radial velocities of the observed sample. The top panel shows 
the distribution for field stars and the bottom four panels show the distributions for cluster 
stars. 
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to a vsini value, an empirical width versus vsini relation is determined by cross- correlating 
many artificially "spun- up" spectra with the same slowly rotating stan dard. These a rtificially 



broadened spectra are constructed using the rotational profiles given in I Gray I (Il992l . e = 0.6). 
This procedure is computed for each order, yielding a f sim estimate for each order. For each 
standard used, a vsini value is computed from the median of all order estimates; a median 
is preferred since the vsini measurements are often upper limits. 

Uncertainties in vsini are determined first from the standard deviation of the multiple 
wsim measurements provided by all standards used. For observations consisting of multiple 
exposures, the uncertainty is the average standard deviation of these measurements. The 
vsini values are assumed to be upper limits if less than 10 km s~^ (just over one-half the 
velocity resolution), or if vsini is less than 3 times the uncertainty in vsini. We find that 
all templates of F, G, and K spectral types give consistent vsini values; an exact match in 
spectral type is not critical in determining vsini. 

As noted in the introduction, poor focus plagued some of our measurements, especially 
on the night of 2002 Feb 2. To assess how this may affect inferred properties, the above 
radial velocity and rotational velocity analysis was conducted on the observations of radial 
velocity standard HD 164922 observed at multiple focus settings (Section 2). Fortunately 
the inferred radial velocities for modest focus offsets, typical of changes over a given night, 
agree with the nominal value to ~ 1 km s~^, consistent with the inferred uncertainties. For 
extreme focus offsets, however, the discrepancy can be as large as 3 km s~^. Similarly, for 
modest focus offsets from nominal, the inferred vsini values would increase modestly, but 
since the standard deviation of measurements among the orders correspondingly increased, 
the end result is simply larger upper limits (e.g. < 16 km s~^ for HD 164922, compared with 
< 10 when observed at nominal focus). However, for extreme focus offsets, the spectrsopic 
blurring resulted in a large apparent vsini (e.g. 22.0 ± 3.3 km s^^ for HD 164922). Since 
the majority of our spectra have at most only modest focus errors, no correction for this is 
made. Rotational velocities inferred from observations made on 2002 Feb 2, however, are 
marked with a colon in Table 1 as likely being artificially too large. 

To empirically assess systematic errors in our radial and rotational velo c ity rn easure- 



ments, our values are compared with values measured by iNordstrom et al. I (120041 . which 
appeared after all of our data were taken). There are 135 stars in common between the sur- 
veys that (1) are not previously identified spectroscopic binaries, (2) are not used as radial 
velocity standards by us, and (3) have radial velocity uncertainties less than 4 km s~^. Figur e 



[H] shows the difference between our radial velocities and those of INordstrom et al. I (120041 ). 
The agreement is quite good; the average difference is -0.2 km s~^ with a standard deviation 
of 2.4 km s~^. With the exception of 2 stars, HIP 67904 and HIP 61072, all measured radial 
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velocity values agree to within 10 km s~^. HIP 67904 and HIP 61072, on the other hand, 
have significantly different radial velocity measurements (> 25cr difference). As discussed 
below, we suggest these stars are spectroscopic binaries. This comparison also demonstrates 
that the radial velocities of stars with vsini values as large as 45 km s~^ can be measured 
without significantly degraded precision. 



In Figure [7] are shown iNordstrom et al. I (120041 ) vs ini values versus the P60 vsini values 
for the overlapping sub-sample described above. While INordstrom et al. I (120041 ) report vsini 
values as low as km s~^ (which is likely unrealistic), our lower resolution data restrict us to 
vsini upper limits of 10 km s~^ or g reater. The maiority of o ur vsini upper limits are thus 
greater than the values measured by iNordstrom et al. I (120041 ). which inhibits interpretation 
of an overall comparison. However, a comparison of the 22 stars with measured values, 
however, sug gests a systematic differe nce. The P60 vsini values are larger, in the median, 
than those of iNordstrom et al. I (120041 ) by 8.6 km s^^. A portion of this difference may be a 
consequence of the spectroscopic blurring caused by the focus errors discussed above. When 
the 3 measured values from 2002 Feb 2 are removed, the P60 values are then larger by only 
6.3 km s~^, in the med .ian. However, an indepe nden t comparison of the 389 st ars with vsini 

Nordstrom et al. (2004) and Strassmeier et al. ( 200oh. upper limits 



measurements in both 
excluded, indicate the 



Nordstrom et al. 



Strassmeier et al. 



(I2OO4J ) values are less than the 
((2000) values by 2.6 km s~^. Thus, a portion of the d ifference with our r aeasu rements may 
be a consequence of artificially low values reported by iNordstrom et al. i (120041 ) . Overall we 
conclude that the majority of our measured vsini values suffer from no systematic biases 
greater than ~ 5 km s~^. The exceptions are stars observed on 2002 Feb 2, which may be 
biased toward artificially large values by ~ 10 — 20 km s~^ because of especially poor focus. 
Since we are uncertain what the systematic correction to the vsini values for these stars 
should be, if any, we make no correction for this in the values presented in Table 1. 



4.3. Temperature and Gravity Indicators 



Within the broad wavelength coverage provided by our echelle spectra are many temperature- 
and gravity-sensitive lines. Empirically, we find that the EWs of the photospheric features 



Fel A4063, Fel A4071, and Call A4226 (see iJaschek fc Jaschek Iil987l ) are especially useful 



temperature diagnostics. In our data, these features all increase in strength in a roughly 
linear fashion from spectral types F8 to K5, before degenerate fiattening and turnover oc- 
curs towards later spectral types. However, fits of the EWs versus B — V color for main 
sequence stars (not shown here) are accurate to only 0.05-0.1 in B — V. This implies that if 
we were to attempt spectral classification from the EWs, spectral type uncertainties of more 
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Fig. 6 . — Difference between P60 radial velocities and values reported by iNordstrom et al 
(l2004j ). The solid squares are likely spectroscopic binaries. 
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Fig. 7. — Comparison of measured vsini values with values from iNordstrom et al. (12004 ) . 

report 



Triangles are upper limits and squares are measured values; INordstrom et al. 



kooi \ 



no upper limits. The 3 x's indicate measurements from the night of 2002 Feb 2, which had 
significant focus problems (Section 4.2). 
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Fig. 8. — Gravity sensitive equivalent width ratios versus B-V color. Dwarf stars are indi- 
cated with solid circles and giant stars are indicated with crosses. 
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than one spectral subclass would result. Therefore we do not use these features to estimate 
spectral types. Spectral types instead are taken from several catalogues of stellar spectral 
classification, when available. These, in order of preference are Houk & Swift (1999), Skiff 
(2005), Kharchenko (2001), Kharchenko, Piskunov, Scholz (2004). The spectral types are 
listed in Table 1; 349 of the 370 stars observed have stellar spectral types. 

For gravit y indicators, we considered the ratios of the gravity sensitive SrII A4216 and 



A4077 lines (see lJaschek fc Jaschek 1119871 ) to the above temperature sensitive lines. In Figure 
IS] we show for our dwarf and giant standard stars the ratios (SrII A4216 + A4077) to (Cal 
A4226 + Fel A4063 + Fel A4071) and SrII A4216 to Cal A4226, as a function of B - V 
color. These ratios are sensitive to both temperature and gravity; the dwarf and giant stars 
are reasonably well separated in these plots. Parabolic fits to the dwarf and giant points 
are shown. The majority of our target objects cluster around the dwarf locus, though some 
are securely between the dwarf and giant loci. Other diagnostic figures, such as Sr II A4077 
/ (Fel A4063 + Fel A4071) versus Cal A4226, (which has the advantage of considering 
only our spectroscopic data), have similar appearances. In principle one might use such 
figures to quantitatively estimate the surface gravities of our program stars; in practice, 
however, the data do not allow this to be carried out with the precision normally expected of 
surface gravity measurements. Rather than estimate surface gravities, we simply identify all 
stars which are more than the la above the dwarf loci in both diagnostic planes. Sixty-six 
potentially low surface gravity stars are identified by these criteria and marked in Table 1. 



4.4. Chromospheric Ca II H and K Activity Measurements 



Core emission in the Calcium H and K spectral lines serves as an indic ator of stellar chro 



mosp heric activity, which is known to be correlated with stellar age (e.g. ISoderblom et al 



19911 ). primarily as a consequence of rotation. We parameteri ze the core emi s sion in our spec- 



tra u sing the standard Mount Wilson Project procedures (IWilson I Il968l : IVaughan et al 



1978[). Our pr e script ion for measuring the core emission closely follows that outlined by 
Duncan et al. I (Il99ll ). First we define the index of core emission, S, as 

where Ni is the number of counts in band i. The continuum bands, R and V, are 20 A wide 
and centered at 4001.07 A and 3901.07 A, respectively. The emission bands, H and K, are 
centered on the emission features whose natural centers are at 3968.470 A and 3933.664 A, 
respectively. The bandpass for the H and K channels is triangular, with a full width at 
half maximum of 1.09 A. Before measuring the counts within these four bands, the spectra 
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are normalized to account for the blaze function and shifted to correct for radial velocity. 
The spectra are normalized in a two step iterative process. First, a ~ 4th order polynomial 
is fit to the spectra; the broad H and K absorption features are excluded in this fit. A 
second polynomial is then fit to the spectra, excluding all spectral features that extend 0.5 
a below the initial polynomial fit. The spectra are then normalized by dividing by this 
polynomial. The velocity shift is determined by cross-correlating the 2 spectra orders with 
the Ca II H and K features, with the spectrum of a bright comparison standard. The implied 
radial velocity is compared with that determined in Section 4.2 to ensure consistency; for 
the chromospheric activity analysis we adopted values determined directly from these orders 
as they typically provide the most accurate shift for these orders. 

With velocity shifted and normalized spectra, the counts within the continuum and 
emission bands of the S index are summed. For the emission bands, the number of counts in 
each pixel is first multiplied by a number from to 1 to mimic a triangular bandpass with a 
full-width at half-maximum of 1.09 A. With these sums, equation ([1]) is used to compute S. 

To correct for systematic offsets in the measured S values relative to previous work, pos- 
sibly due to slight diffe rences in the method ology or data quality, the S values are compared 



to those measured by iWright et al. I (I2004J . also appearing after all data for this program 
were gathered). These comparisons permit transformation of the S values to what we refer 
to as Mount Wilson S values, Smw- This is realized as a first order regression analysis on 
log 5* in the form of 

log Smw = A + B\ogS. (2) 
The uncertainty, cfsmwi given by the standard formula 

2 fdSMW V I f'^^MW V to\ 

^s..= [-^-a) +[-^-b) . (3) 

The analysis was performed on 8 of the 9 observing runs and for the entire overlap- 
ping sample; there were too few overlapping stars for the 2001 December observing run for 
comparison. Fortunately, the uncertainties in the implied correction estimates for 7 of the 8 
observing runs agree well with that determined for the entire overlapping sample. The one 
exception was for the July 2001 observing run, for which the corrections appeared signifi- 
cantly different. Therefore, the corrections determined for this run were used to transform 
its 5" values to Smw values, while the remaining runs used the corrections determined for the 
entire sample, which are more precise. The resulting Mount Wilson Calcium index values 
and uncertainties are listed in Table 1. Stars with uncertainties that are possibly much larger 
than reported, because of a large systematic uncertainty that is not formally propagated (e.g. 
radial velocity error) are marked with an "e" in Table 1. 
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Finally, the Smw values are converted to R'hk values, which express the activity as a 



fractio nal ratio with the underlying star, following the prescription described in lNoyes et al 
( ll984j ). However, this conversion requires proper characterization of the underlying photo- 
spheric flux, which can be assessed from the intrinsic stellar B — V co lor. B — V colors 



for most stars are obtained from the Tycho-2 catalogue (lH0g et al. Il2000f). These valu es are 



transformed to Johnson B — V colors, following the relations in lMamajek et al. I ( 2002 1^. For 
objects without Tycho photometry, the B — V colors are est imated from the stellar spectral 
type (Section 4.3) following the relations of I Johnson I (119661 ). 



R'hk is calculated as the difference in the logarithms of the flux calibrated core emission, 
Rhk, and photospheric emission, Rphot'- 



R!. 



HK 



Rhk — R 



■phot- 



(4) 



Here Rhk is calculated from the Smw values and calibrated using the stellar B — V color, 
and Rphot is calculated from the B — V color alone (see Appendix b in Noyes et al. 1984 for 
explicit formulae). Although the calibration of these terms is determined over the restricted 
B — V range of 0.44 < {B — V) < 0.82, we use it for the full B — V range of the observed 
sample (0.19 < {B — V) < 1.65). The resulting R'^k values are hsted in Table 1. 



5. Sample Properties 

5.1. Spectroscopic Binaries 

Table 2 lists the 8 stars with photospheric spectral features that are doubled in our 
observations. Two of these double-lined spectroscopic binaries were previously known (HE 
848, SAO 178 272), while HD 78899, TYC 7310 503 1, TYC 7305 380 1, and HD 140374 are 
newly identified. The multiple epoch observations of some stars also allows us to identify 
single-lined binaries based on changes in radial velocity. We consider an object a candidate 
binary if its radial velocity changed by more than lOcr between observations; 3 stars meet this 
criteria: HD92855 (HIP52498), HD132173 (HIP73269), HD61994 (HIP38018; see T a ble 1) . 



Similarly, based on comparisons with radial velocities published in lNordstrom et al. I (120041 ). 
HD108944 (HIP61072) and HD121320 (HIP67904) are candidate binaries (see FigureE]). All 
5 are newly identified candidate binary systems. Finally, since the radial velocity dispersion 
in young clusters is typically < 2 km s~^, comparison of cluster member radial velocity 



^There is a typographical error in the relations of Mamaiek et al. ( 20021 ). The +0.007813 in equation 
C6 should be -0.007813. 
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measurements with the median cluster values can identify single-lined spectroscopic binaries 
Based on these comparisons, the two Pleiades members HIT 571 and LH9 8 103 are identified 



as spectroscopic binaries; both have been identified previously as such (IMermilliod et al. 



19921 : iQueloz et al. Ill998f ) 



5.2. Sample Age Dispersion 

Several of the quantities extracted from the observed spectra are likely to be correlated 
with stellar age, based on extensive studies in the literature. Generally speaking, younger 
stars tend to be more rapidly rotating, have larger lithium abundances, and have stronger 
chromospheric emission lines such as Call H&K and Ha. Thus the measured quantities 
vsini, EW[LiI A6708], log Rhk^ EW[Ha] should correlate with age. Unfortunately the 
precision with which stellar ages can be estimated from these quantities is in general poor, in 
part because of the additional dependencies upon stellar spectral type and rotation. Rather 
than attempting to disentangle these complicating effects for individual stars, we simply 
present the data to illustrate the empirical correlations amongst the measured quantities 
and to assess the overall age range of the observed sample. 

Figure shows the distributions of EW[LiI 6708], EW[Hq;] and t>sim versus logi?^^. 
Error bars are not displayed due to crowding, but can be inferred from the information in 
Table 1. The EW[LiI] and log R'hk values are reasonably well correlated, as expected if 
strong lithium absorption and strong chromospheric emission, as traced by Ca II emission, 
both indicate stellar youth. The EW[Ha] and log R'^k values are at best weakly correlated; 
stars with larger log R'^^^ values have, on average, smaller EW[Ha] values, suggesting they 
are partially filled in by emission. The strong spectral type dependence of this temperature 
sensitive line likely masks much of the underlying correlation which might be revealed by 
considering an analogously defined quantity R'^a- The vsini values are well correlated, in 
the mean, with log R'jj^ values. This is expected since increased rotation is believed to cause 
increased chromospheric activity, a characteristic common in young stars. Finally, in Figure 
[To] are shown the distributions of EW[Ha] and vsini versus EW[LiI]. Both the EW[Ha] and 
fsim values are correlated with EW[LiI], in the mean, though both show a large scatter at 
all EW[LiI] values. 

The above correlations corroborate the suggestion that these measured quantities trace 
stellar age. While the large scatter in these correlations inhibits assigning specific ages to 
individual stars, the ensemble distributions can help assess the age of the observed sample. 
Since the stellar temperature can also affect these quantities, we consider these distributions 
as a function of B — V color, a proxy for stellar temperature. The range of observed values at 




Fig. 9. — EW[LiI 6708], EW[Hq;] and vsini versus the chromospheric activity diagnostic 
logi?^^. Measured values are shown as squares while upper limits are shown as triangles. 
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Fig. 10. — EW[Hq;] and vsini versus the EW[LiI 6708]. Measured values are shown as 
squares while upper limits are shown as triangles. 
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any B — V color should therefore (mostly) be a consequence of the age spread in the sample; 
a comparison of the spread of values at different colors can serve as a check for age biases 
with respect to temperature. 

In Figure [TT] are shown the distributions of t>sim, EW[Li I A6708], log R'hk^ EW[Ha] 
versus B — V color (as assembled in Section 4.4). The observed dwarf standards are distin- 
guished in order to illustrate the properties of old slowly rotating stars. The distributions of 
vsini and log R'fjx are independent of B — V color, indicating a similar spread in rotation 
rate and chromospheric activity over the range of spectral types observed. The EW[LiI] 
values span a broad range at each B — V color. The larger EW[Li I] values at the red end is a 
temperature effect; f or a specified lithium abunda nce, the EW[LiI] increases with decreasing 



stellar temperature (jPavlenko fc Magazzu Ill996l ). The EW[Ha] values show the strongest 



B — V trend, as expected for this temperature sensitive line; the EW[Ha] values decrease 
with increasing B — V color. There is also a larger scatter at larger B — V values. The 
scatter above the dwarf loci can be interpreted as stars with some Ha emission, whether 
this is partially or completely (e.g. EW[Ha] < 0.0) filling in the photospheric absorption 
line. Four stars with strong Ha emission (EW[Ha] > —10 A) are not plotted; their profiles 
are shown in Figure |H Stars with some Ha emission span the full range of B — V colors, 
though there appear to be many more at redder colors [B — V > 1.0) than bluer colors 
[B — V < 0.6). However, since an EW is defined relative to the continuum flux, a unique 
EW[Ha] value corresponds to much less Ha flux for cooler (red) stars than hotter (blue) 
stars. Overall, these age indicators are consistent with the observed sample having ages 
that span from a few million years to a few billion years, as was intended. Moreover, all 
spectral types appear to span this range of ages; there is no strong age bias with spectral 
type amongst the observed sample. 



5.3. Call H&K Saturation in Young Stars 



While log R'fjK values have been shown to correlate with stellar age (e.g. lSoderblom et al. 



199ll ). this diagnostic typically has not been used to estimate ages less than ~ 100 Myr. This 
primarily stems from the lack of R'hk measurements for stars in young clusters of known 
age, from which the age relation could be calibrated. Although the P60 R'hk measurements 
of young cluster members presented here may help establish this calibration, we illustrate 
that accurate ages will be inhibited by the broad range in rotation rates for stars at a given 
young age and the saturation of flux at high rotation rates. 

In Figure [12] are shown the R'^k measurements versus projected rotational velocity 
(fsini) of stars in the open clusters represented in our sample, including a Persei, the 
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Fig. 11. — Projected rotational velocity (vsini), EW[Li A6708], log R'jjx EW[Ha] versus 
B — V color. Filled circles indicate dwarf standard stars. For clarity, a few stars in the 
observed sample are bluer than (SAO 41857) or redder than (HIP 108467, HIP 90190, and 
HIP 60783 ) the B-V range shown here, the values are nevertheless included in Table 1. 
Similarly, in the bottom panel four stars with strong H« emission (EW[H«] < —10; Figure 
4) are not shown. 
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Pleiades, the Hyades and collectively (post) T Tauri stars in Upper Scorpius, Upper Cen- 
taurus Lupus, and R Corona Australis. As expected, the R'hk measurements are correlated 
with f sinz; the most rapidly rotating stars have, on average, the largest R'hj<;; values. How- 
ever, Figure [12] also indicates that once the fsim values increase above ~ 30 km s~^, the 
R'hk values stay roughly constant, implying no increase in the Call H&K emission line flux 
toward larger rotation rates. A similar "saturation" phenomenon has been seen before based 
on observations of other emission lines originating in ste llar chromos pheres (01 1304, CII 



1335), as well as X-rays originating in stellar corona (e.g. IVilhu Ill984l ). It is not yet known 



whether the cause for flux saturatio n is internal (e.g. dynamo) or external (e.g. active re- 



gions) to the star (see iGiidel II2004I ). Saturation of the Call H&K fluxes, as measured by 



the R'hk diagnostic, has not been demonstrated before. 

While saturation of the Call H&K flux seems consistent with other chromospheric 
and coronal observations of rapidly rotating stars, we considered the possibility that this 
saturation is simply a consequence of the method used to determine the S values, which are 
used to calculate R'hk- Speciflcally, since the prescription outlined in Section 4.4 involves 
measuring the Ca II K and H emission with a triangular passband of specifled width, apparent 
saturation could be caused by emission line flux being rotationally broadened outside of this 
passband, and thus not included in the measurement. To test this, we remeasured S values 
for a set of stars that were artiflcially rotationally broadened from 10 to 150 km s~^; we 
assume the core emission is rotationally broadened by the same amount as the photosphere. 
This showed that the technique used to measure the S values begins to miss core emission 
flux above vsim of ~ 100 km s~^. Since the saturation observed for R'^^ occurs well below 
this value, we interpret the saturation as a real effect. 

To further demonstrate the case for saturation, in Figure [T2] we also illustrate the ratio of 
X-ray flux to bolometric flux (Lx/Lhoi) for the stars in these clusters which have the available 
data. X-ray fluxes are from ROSAT observations and bolometric fluxes are determined from 
the atmospheric model which best flts the stellar energy distribution (see Mamajek et al. 
2007 for a complete description of the Lx / derivation). A saturation trend is identifled for 
coronal X-rays similar to that for chromospheric Ca II H&K core emission. The most rapidly 
rotating stars are the most fractionally X-ray bright, but Lx/Li,oi remains roughly constant 
above vsini values of ~ 30 km s~^. This phenomenon has been demonstrated previously by 
many authors. What we have newly demonstrated here is that the same young cluster stars 
that have saturated X-ray emission also have saturated chromospheric Ca II H&K emission. 



We are grateful to the observing assistants formerly present at the Palomar 60-inch, 
Karl Dunscombe and Skip Staples, for their help and patience in acquiring this large data 
set. We thank E. Mamajek and C. Slesnick for assisting with the observations and Jonathan 
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Fig. 12. — R'hk ^"^^ Lx/LboI measurements of open cluster members versus vsmi. Filled 
circles are a Per members, squares are Pleiads, triangles are Hyads and open circles are T 
Tauri stars in Upper Scorpius, Upper Centaurus Lupus, and R Corona Australis. Stars with 
a vsmi upper limit are plotted with a vsim of 5 km s~^ for clarity. 
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Table 1. Spectroscopic Properties of Observed Single(- Lined) Stars 



HD 


HIP 


other 


Grp 


RA(2000) 


DEC (2000) 


Epoch 


SNR 


EW[Li I] 


EW[Ha] 


RV 


vsini 


SNR 




Spectral 
























atLi 


(A) 


(A) 


(km 


(km 


atHK 




Type'' 


B-V'' 


224873 


110 


RX ,10001.4+3936 


Field 


00 


01 


23.66 


+39 36 38.1 


2001 


Dec 2 


105 


0.032 


1.74 


-5.6 ± 1.4 


< 16.3 


50 


0.365 0.068 


KO 


0.76t 


106 


490 


RX J0005. 9-4145 


Field 


00 


05 


52.55 


-41 45 10.9 


2001 


Dec 2 


217 


0.156 


3.06 


1.6 ± 1.2 


< 18.2 


16 


0.675e 0.118 


GOV 


0.59t 


377 


682 


IRXS J000826.8+063712 


Field 


00 


08 


26.74 


+06 37 00.6 


2001 


Jul 27 


348 


0.122 


2.34 


-3.8 ± 2.2 


< 12.0 


135 


0.348 0.12 


G2V 


0.63t 


691 


919 


V344 And 


Field 


00 




22.40 


+30 26 68.4 




Oct 29 


266 


0.122 


2.20 


-2.4 ± 0.8 


< 10. 


84 


0.388 0.072 


KOV 


0.78t 


















00 


Sep 18 


174 


0.115 


2.14 


-3.0 ± 0.8 


< 13.7 


84 


0.391 0.073 




0.78t 


984 


1134 


IRXS J001410.1-071200 


Field 


00 


14 


10.25 


-07 11 57.0 


2001 


.Tul 27 


315 


0.099 


2.99 


-3.2 ± 2.2 


42.1 ± 1.7 


131 


0.296 0.104 


F7V 


0.50t 


1326A 


1475 


G.7 15 A 


Mstar 


00 


18 


22.57 


+44 01 22.2 


2002 


Oct 30 


198 


< 0.01 


0.51 


12.0 ± 0.9 


< 20.5 


38 


0.568 0.101 


M2V 


1.48t 


132fiB 




GJ 15 B 


Mstar 


00 


18 


25.50 


+44 01 37.6 


2002 


Oct 30 


87 


< 0.02 


0.17 


10.7 ± 0.9 


< 18.9 


14 


1.083 0.181 


M5V 


1.61s 


3765 


3206 


GJ 28 


Field 


00 


40 


49.29 


+40 11 13.3 


2002 


Sep 17 


259 


0.015 


1.31 


— 63.8 ± 0.7 


16.3 ± 5.1 


65 


0.218 0.043 


K2V 


0.97t 






QT And 


Field 


00 


41 


17.32 


+34 25 16.8 


2002 


Oct 29 


160 


0.128 


-0.17 


6.0 ± 0.8 


24.4 ± 2.2 


40 


1.617 0.259 


G 


0.94t 


















2002 


Sep 18 


83 


0.166 


-0.23 


5.1 ± 0.7 


26.1 ± 2.2 


27 


0.705 0.123 




0.94t 


6434 


5054 


GJ 9037 


Field 


01 


04 


40.15 


-39 29 17.4 


2001 


Dec 2 


231 


< 0.004 


2.37 


21.9 ± 1.2 


< 10. 


28 


0.396 0.074 


G2/3V 


0.60t 


6963 


5521 


SAO 36995 


Field 


01 


10 


41.92 


+42 55 54.7 


2002 


Sep 17 


206 


< 0.003 


2.05 


-31.7 ± 0.7 


< 12.58 


77 


0.226 0.045 


G7V 


0.73t 


7661 


5938 


SAO 147702 


Field 


01 


16 


24.18 


-12 05 49.3 


2003 


Feb 8 


185 


0.075 


2.04 


5.2 ± 0.4 


< 10. 


33 


0.479 0.087 


KOV 


0.77t 




6276 


SAO 147747 


Field 


01 


20 


32.26 


-11 28 03.6 


2003 


Feb 8 


178 


0.153 


1.99 


8.3 ± 0.4 


< 10. 


29 


0.619 0.110 


GO 


0.80t 


8467 


6676 


SAO 64666 


Field 


01 


24 


27.98 


+39 03 43.6 


2002 


Sep 17 


172 


< 0.011 


1.64 


15.0 ± 0.7 


< 10. 


60 


0.278 0.054 


G6 


0.77t 


8941 


6869 


SAO 92463 


Field 


01 


28 


24.36 


+17 04 46.2 


2001 


Jul 29 


421 


0.006 


2.59 


9.3 ± 2.2 


< 10.96 


181 


0.175 0.035 


F8IV-V 


0.62t 


















2002 


Sep 18 


233 


< 0.004 


2.69 


8.29 ± 0.70 


< 12.3 


123 


0.190 0.069 


F8IV-V 


0.62t 


8907 


6878 


IRAS 01256+4200 


Field 


01 


28 


34.35 


+42 16 03.8 


2001 


Jul 26 


266 


0.048 


1.55 


6.3 ± 2.2 


< 14.34 


142 


0.263 0.093 


F8 


0.49t 


9472 


7244 


SAO 74789 


Field 


01 


33 


19.03 


+23 58 32.1 


2002 


Sop 18 


190 


0.055 


2.28 


10.9 ± 0.7 


19.33 ± 4.22 


83 


0.277 0.053 


GO 
23 K3Ve 
KOV 


0.68t 






RX J0137.6+1835 


Field 


01 


37 


39.41 


+ 18 35 33.2 


2001 


Jul 28 


194 


0.416 


-0.59 


1.8 ± 2.2 


21.25 ± 3.26 


26 


0.512e 0.171 


1.03t 


10780 


8362 


IRAS 01441+6336 


Field 


01 


47 


44.88 


+63 51 09.1 


2001 


Jul 28 


203 


0.012 


1.66 


1.4 ± 2.2 


< 10. 


66 


0.237 0.085 


0.80t 


11860 


9073 


SAO 76038 


Field 


01 


66 


47.28 


+23 03 04.1 


2002 


Oct 29 


176 


0.018 


2.32 


3.3 ± 0.8 


< 12.08 


64 


0.287 0.055 


1 G6 


0.69t 


















2002 


Sep 18 


138 


0.016 


2.27 


1.6 ± 0.7 


< 12.38 


76 


0.292 0.056 




0.69t 



^Stars with especially uncertain S values are marked with an "e". 
^Stars which appear low surface gravity are marked with an asterisk. 

*^B-V colors from the Tycho catalogue are marked with a "t" while those from the stellar spectral type are marked with an "s". 
*^Stars with B-V colors outside the range over which R'hK" values are calibrated are marked with a "c". 
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Table 2. Observed Double-Lined Binary Stars 



SNR EW[Lii] EW[C;aI] BW[Ha 
HD HIP other Assoc RA(2000) DEC(2000) Epoch atLi (mA) (mA) (A) 







CI Molotte 20 848 


aPer 


03 


29 


26.24 


+48 12 11.74 


2002 


Oct 31 


149 


45/69 


61/85 


1.65/0.95 


78899 


45187 


GJ 3539 


Field 


09 


12 


28.24 


+49 12 23.4 


2003 


Feb 9 


260 


11/50 


91/121 


0.83/0.72 






SAO 178272 


Field 


09 


59 


08.42 


-22 39 34.57 


2002 


Feb 4 


107 


87/162 


105/153 


em. 


















2001 


Dec 2 


134 






em. 


















2002 


Apr 19 


51 


103/204 


174/233 


-0.1/-0.2 






TYC 7310 503 1 


Field 


14 58 


37.69 


-36 40 30.27 


2003 


Jun 5 


82 


121/222 


99/180 


em. 






TYC 7305 380 1 


Field 


14 


50 


25.82 


-35 06 48.66 


2002 


Feb 3 


85 


132/250 


103/66 




















2003 


.Tun 5 


79 


307/60 


120/54 




140374 


77081 


TYC 7331 1235 1 


Field 


15 


44 


21.06 


-33 18 54.97 


2003 


,Jun 4 


145 


112/116 


61/77 


-0.96/-1.17 


142229 


77810 


SAO 121238 


Field 


15 


53 


20.01 


+04 15 11.5 


2002 


Feb 2 


119 


56/121 


29/79 




218738 


114379 


KZ And 


Field 


23 


09 


57.34 


+47 57 30.0 


2003 


Jun 6 


73 


69/107 


108/165 





